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Abstract 

The soliton formation is considered in the Nambu-Jona-Lasinio model with 
local four quark interaction and various schemes to regularize the energy con- 
tribution of the polarized vacuum. No additional constraints are admitted in 
order to stabilize the soliton. While solitons are unstable in the proper-time 
regularized version the three momentum cut-off regularization apparently is 
more appropriate. Using a semi-classical approach multi-quark solitons ob- 
tained from that scheme are discussed. However, no self-consistent non-trivial 
unit baryon number configuration has been found. We also study a renormaliz- 
able extension of the model. In this case no stable multi-quark solitons are 
obtained within the semi-classical approach. 
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1 Introduction 



The strong interaction dynamics of mesons and baryons is described by quantum 
chromodynamics (QCD). It is assumed that this theory explains hadronic phenom- 
ena in terms of the quark and gluon degrees of freedom. At low energies, comparable 
with the low lying hadron masses, QCD exhibits a non-perturbative behavior. This 
circumstance renders the analytic study of the theory rather difficult. Nevertheless, 
a qualitative description of important aspects of QCD is possible on the basis of 
effective theories which have in common with QCD the symmetries of the quark- 
flavor dynamics, in particular the pattern of spontaneous chiral symmetry breaking. 
In this respect the Nambu-Jona-Lasinio (NJL) model has been very successful 
in the description of the meson sector 0, [J. 

It is therefore natural to investigate the soliton solutions of the model in order 
to describe baryons according to the 1/Nc expansion. In this context the important 
issue to be addressed concerns the stability of these soliton configurations. To be 
specific we consider the NJL model with a local four quark interaction and various 
regularization prescriptions for the contribution of the polarized vacuum to the mass 
of the soliton. This is motivated by the common belief that its stability strongly 
depends on the imposed regularization procedures. For example, in refs ||, || the 
description of baryons as solitons with three valence quarks derived from the NJL 
model with proper-time regularization || is reviewed. For this description to work 
the scalar and pseudoscalar fields a and tt must artificially be stabilized by the so- 
called chiral circle condition a 2 + 7r 2 = const.. However, it is known that without this 
condition the soliton using proper-time regularization is not stable @, §]. Here, we 
wish to define the stable soliton as a non-trivial mean-field solution, which is stable 
by itself once the regularization prescription has been fixed. In particular we do not 
allow for additional constraints which artificially stabilize the field configuration but 
do not directly follow from the model. This criterion is apparently not met by the 
just mentioned soliton. 

Another example is the application of the proper-time regularization to the 
imaginary (anomalous) part of the effective action as well 00. As a consequence 
the baryon number is no longer quantized and may assume non-integer values. The 
soliton can then be stabilized by constraining the regularized baryon number e.g. to 
unity. As this stabilization mechanism originates from a particular regularization 
prescription the question of whether or not it fulfills our above criterion remains a 
matter of taste. Alternatively one might add further interaction terms to the model 
Lagrangian. In that respect the inclusion of a Higgs-type meson self-interaction, 
i.e. (a 2 + 7r 2 — const.) 2 , has been shown to render the soliton stable [TO. 

Not to be misinterpreted, we do not regard these additional constraints as un- 

3 This term is conditionally finite and hence does not require regularization. Its regularization 
might spoil the anomaly structure of the model. 
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reasonable. The chiral circle condition as suggested by Skyrme type models, which 
have this condition implemented by definition, may certainly be in accordance with 
a reasonable description of single baryons as solitons. The baryon number constraint 
to integer values is obvious, and the forth-order self-interaction is unavoidable in 
sensibly renormalizable models. An additional forth-order self-interaction may also 
be motivated by higher multi-quark interactions following from QCD when gluons 
are integrated out [fUH - The point, however, is that these additional constraints 
do not directly follow from the NJL model together with an adopted regularization 
procedure. 

Recently it has been shown [fl^| , using the Thomas-Fermi method, that the 



NJL model with three dimensional momentum cut-off regularization of the Dirac 
sea, apparently accounts for the existence of stable solitons, provided the number 
of valence quarks is large enough. Needless to emphasize that the phenomenology 
of multi-quark or multi-baryon systems requires the a and 7r fields to be treated 
independently. We compute the critical size of this soliton below which it ceases to 
be stable. This size is determined by the interplay between the bulk and surface 
energy densities. For the proper description of the surface the lowest order gradient 
corrections to the Thomas-Fermi method in the Wigner-Kirkwood expansion are 
taken into account. These investigations are described in section 2 and are completed 
in section 3 with a self-consistent calculation for unit baryon number imposing the 
corresponding regularization scheme on the single particle Dirac Hamiltonian. 

Finally, in section 4, we briefly discuss a renormalizable extended version of the 
NJL model whose solitons are yet unstable, apparently for the same reason as those 
in the proper-time regularized version. Our conclusions are summarized in section 5. 



2 Thomas— Fermi method and Wigner-Kirkwood expansion 

Our starting point is the semi-bosonized version of the NJL model where the 
scalar and pseudoscalar fields a and tz have been introduced through a Hubbard- 
Stratonovich transformation to eliminate the original four-fermion interaction 

£njl = q [i'fdft - (a + ij 5 r -n)]q . (1) 

We do not consider explicit chiral symmetry breaking. In the vacuum sector the 
scalar field is identical to the constituent quark mass m. In order to prevent the 
model from becoming a trivial theory of non-interacting mesons a regulator has to 
be retained or alternatively the necessary counterterms have to be added to the La- 
grangian ([]]). In the subsequent section we will study the three-momentum cut-off 
in more detail because this regulator suppresses the high momentum components, 
which in the proper-time regularization scheme cause the soliton to shrink to a 
point-like singularity |7], |j. For that reason we present the following expressions in 
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that particular scheme, the generalization to other regularization schemes is straight- 
forward. The three-momentum cut-off A is related to the four-fermion coupling G 
and the constituent quark mass m via the gap-equation 

d 3 p 1 



1 

G 



AN, 
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(27r) 3 y/p' 2 + m 2 
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In the chiral limit the pion decay constant obtained from 

d 3 p 1 
(2vr) 3 (p 2 + m 2 ) 3 / 2 



fl = N c m 2 J 



(3) 
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will later be used to determine the value of the constituent mass m as a function of 
the cut-off A. 

The semi-classical, relativistic particle and energy densities for fermions in the 
background of a scalar field are derived within the Wigner-Kirkwood (WK) expan- 
sion up to order h 2 0, 0, || 
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The leading order (Thomas-Fermi method) was already investigated in Ref . [I2j and 
the next to leading order includes gradient corrections 
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Note that the scalar field is space-dependent but does not depend on the momenta. 
For the valence contribution the three-momenta were integrated from zero to the 
Fermi momentum pp. To this the contribution of the polarized vacuum has been 
added. It is obtained from the same integrals, however, with the opposite sign and 
the upper bound changed to the cut-off A. This not only explains the integra- 
tion bounds in eq. (^) but also removes the infra-red singularities dwelling in the 
individual contributions to the gradient terms. 
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In principle pseudoscalar fields may be introduced in a chirally symmetric way by 
the replacements a 2 — > <7 2 + 7r 2 , (Vcr) 2 — > (Vcr) 2 + (V"7r) 2 and crAcr — > crAo- + 7r • A7r. 
We do not give the explicit expressions because within the semi-classical approach 
the pseudoscalar fields are not excited. The densities (H) 

p = p- p( mc ) , 



£ = £- £ {vac) 

determine the quark number and the energy measured relative to the vacuum. All 
further applications follow from the variation of the functional 

5 J d 3 r{£ - Xp) = , (7) 

where the Lagrange multiplier A (chemical potential) has been introduced to fix the 
quark number A. This constraint should not be confused with the one employed in 
[[]]. Here it is due to the semi-classical treatment and serves solely to determine the 
Fermi momentum pp. In a microscopical calculation as e.g. the one presented in 
subsection 3.3 no such constraint is required. 

Subsequently we apply this method to hadronic matter and to finite hadronic 
systems. 



IpI<A V 

m 2 f d 3 p 



IpI<A 



3 Three dimensional momentum cut— off regular izat ion 

Apparently the three-momentum cut-off regularization for the vacuum contribu- 
tions is a good candidate to yield stable solitons. It successfully suppresses the high 
momentum components contained in the infinitely high and narrow peak of the 
scalar field at the origin which arises e.g. in the proper-time regularization scheme 
[0, HJ. In the following subsections we apply the Thomas-Fermi method together 
with its gradient corrections to hadronic matter and to finite hadronic systems. Fi- 
nally, in subsection 3.3, we present a fully self-consistent calculation for unit baryon 
number. 



3.1 Hadronic matter 

For hadronic matter the gradient terms in (|j) do not contribute, only the leading 
order in the WK expansion which corresponds to the Thomas-Fermi method sur- 
vives. The variation (|7|) with respect to the Fermi momentum pp and the scalar 
field a, respectively, leads then to simple expressions for the particle number and the 
energy discussed already in [12] in connection with the chiral phase transition. The 
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binding energy per quark S/p — m is plotted in Fig.l of that reference for various 
ratios of the cut-off in units of the constituent mass. The maximal binding energy 
(minima of those curves) is always reached in the symmetric phase, a = 0, at 



£/p — m=p F — m 
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which in units of the free constituent quark mass m is the function of A/m depicted 



Figure 1: Binding energy per quark for hadronic matter in dependence of the cut-off 
A. Solitons are obtained for A/m < 1.80. 

in our Fig.l. From this figure we notice that hadronic matter becomes bound for 
A/m < 1.80, for larger ratios no solitons exist. Since for finite systems the gradient 
terms (cf. next subsection) introduce surface repulsion, the latter statement is quite 
general. With decreasing ratio A/m we expect to find solitons also for a smaller 
number of quarks. For what follows we choose A/m = 1.7, 1.6 and 1.5 which corre- 
sponds to a NJL coupling G = 1.60, 1.86 and 2.18 m 2 according to the gap-equation 
(0). All dimensionalfull quantities will be expressed in terms of the free constituent 
mass m. This scale may be fixed by the empirical value of the pion decay constant 
f n ~ 90MeV in the chiral limit using (Q). The corresponding model parameters and 
the results for the binding energy per quark in hadronic matter are listed in Table 
1. The last row refers to an unphysically low value for the cut-off A/m = 1.0 with 
strong coupling G = 6.17m 2 leading to a considerable increase of the binding energy 
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Table V. Parameters of the model for various ratios of the cut-off divided by 
the constituent mass. The pion decay constant is kept fixed at f n = 90MeV. The 
last two columns give the quark condensate and the binding energy per quark in 
hadronic matter. 



A/m 


m [MeV] 


A [MeV] 


G [GeV" 2 ] 




(- < qq >)^[MeV] 


B/p-m [MeV] 


1.7 


348 


592 


13.2 


298 


-5.8 


1.6 


365 


583 


14.0 


297 


-12.5 


1.5 


383 


575 


14.8 


296 


-20.6 


1.0 


553 


553 


20.2 


301 


-101.8 



per quark. We provided this parameter set for later reference in subsection 3.3, 
where we discuss the three quark soliton as obtained in a self-consistent calculation. 

In the following subsection we are taking the next to leading terms in the WK 
expansion into account in order to study the soliton formation in finite systems. 

3.2 Finite hadronic systems 

Here we consider the second order terms in the WK expansion (|5|), whose gradi- 
ent terms contribute a repulsive mesonic kinetic energy corresponding to the surface 
repulsion of the soliton. Since this semi-classical expansion is designed for large par- 
ticle numbers, spin-isospin uncorrelated quark states are assumed. Consequently the 
pseudoscalar field enters these expressions at least quadratically and the variation 
(0) makes this field vanish identically. Note that the non-zero value of the scalar 
field in the vacuum provides the driving term for cr(r). For symmetry reasons it is 
a radial function. 

The variation (|7|) with respect to the Fermi momentum Pf{t), which is now a 
radial function as well, leads to 

e F = Jpf + c 2 = A r < R 'inside' 

(l-3a 2 /p 2 F )(Va) 2 -2aAa = 2Ap F r>R 'outside'. (9) 

The "radius" R divides coordinate space into an interior region, where the Fermi 
energy is fixed by the chemical potential as in the hadronic matter case, and into an 
outer region, where in the Thomas-Fermi approximation pp = when the gradients 
are neglected. The presence of these gradient terms prevents pp and hence the quark 
density from becoming exactly zero in the outer region. Instead all these quantities 
will obtain a smooth exponential tail. 
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In the interior region a non-linear differential equation follows from variation (|7|) 
with respect to u{r). For completeness we give this equation explicitly 




where we have introduced the abbreviation e\ = VA 2 + a 2 . This second order 
differential equation, subject to appropriate boundary conditions, is then solved for 



a given chemical potential A = ep(R) = yp F (R) + & 2 (R) related to a corresponding 
quark number A. In practice this is achieved by fixing the radius parameter R. 

In principle, a similar differential equation has to be solved in the outside region 
together with @, resulting in an additional set of two coupled non-linear second 
order differential equations for the functions a(r) and pF{r)- Fortunately for large 
distances the asymptotical solution is known analytically 



The constant appearing in cr(r) is determined by the boundary conditions at r = R. 
We assume these solutions to be valid in the entire outer region, which is justified 
by the observation that at the radius R the profile cr(R) is already close to its 
asymptotical value m (e.g. for A/m = 1.5 with A = 12 quarks we have cr(R) = 
0.943m at R = 3.12m -1 , cf. solid curve in Fig.3). Using a logarithmic scale the 
resulting binding energy per quark £/p — m is plotted in Fig. 2 as a function of the 
quark number A for various values of the three-momentum cut-off. The binding 
is always weaker than that of hadronic matter which is slowly approached with 
increasing particle number. The surface effect increases rapidly with decreasing 
particle number such that there exists a critical number of quarks for which the 
soliton ceases to exist. 

This number sensitively depends on the ratio A/m which controls the coupling 
constant G and the strength of binding in the NJL model with three-momentum 
cut-off. This dependence is shown in more detail in Fig.3 where the quark number is 
again plotted logarithmically. For A/m > 1.8, as we have seen above, quark matter 
becomes unbound. With decreasing ratio then also solitons with a smaller number 
of quarks become bound, e.g. for A/m = 1.5 we may expect solitons with A > 10 





■oo 



m 3 (m — cr(r)) , 



(11) 



8 



2 

log A 



Figure 2: Binding energy per quark as a function of the number of quarks A. The 
solid, dashed and dotted lines corresponds to different values of the cut-off A/m = 
1.5, 1.6 and 1.7 respectively. The corresponding infinite matter binding energies 
(Table 1) are -0.054, -0.034 and -0.017m. 



1 .4 
cutoff [m] 



Figure 3: Critical number of quarks in dependence of the cut-off A/m as obtained 
in the semi-classical approximation. Above the critical curve soliton formation is 
possible. The dashes indicate that this result should not be trusted for small quark 
numbers. 
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quarks to exist. However, we should be cautious in trusting our results for too small 
quark numbers. Already for A = 12 with A/m = 1.5 (solid curve in Fig. 2) the 
second order WK contribution to particle number and energy respectively, amounts 
to 20% and we may expect an error of about 5% from neglecting the higher orders 
in this expansion. In particular, for the case of special interest, A = 3, a hedgehog 
solution with non-vanishing pseudoscalar field is required. This should alter the 
results considerably. Finally in Fig. 4 we show the calculated profiles for A/m = 1.5 



0.8 
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0.4 
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1 1 "1 1" 1 1 1 1 1 

01 2345678 

radius [1/m] 

Figure 4: Profiles a for A/m = 1.5 as function of the radius in inverse quark masses. 
The solid, dashed and dotted lines represent A = 12, 60 and 120 quarks. 

for various quark numbers. For large particle numbers (A = 60, 120 etc.) the profiles 
differ only by their radii. Inside they are essentially zero and outside they approach 
the free constituent mass, both regions being smoothly connected at the surface of 
the soliton. For small particle numbers (A = 12, solid line) the attraction is no 
longer strong enough to drive the scalar field into the interior region to zero, instead 
the profile starts in the origin at a finite value. 

3.3 Self— consistent calculation for A = 3 

Here follows a short description of searching for a self-consistent configuration with 
unit baryon number in the three momentum cut-off regularization scheme. For this 
purpose we consider the following energy functional 

£[*, *] = ^vai (1 - sign(e val )) - ^ £ e„(A) + ± J d 3 r (a 2 + tt 2 - m 2 ) . (12) 
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In contrast to other covariant regularization schemes which regularize the functional 
trace, which arises when integrating out the fermions, like e.g. proper-time [|, ||, 
the three momentum cut-off scheme is already applied at an earlier stage, namely 
at the level of the single particle Dirac Hamiltonian which is projected onto the 
corresponding subspace in momentum space. That is, the single particle energies 
e M (A) are eigenvalues of 

h A = P A (a-p + p[a + t l5 T-7z])P A , (13) 

where the projection operator acts on the spinorial wave-functions, 

pm -«j J 1 for (/i|p 2 |/i) < A 2 , , 

Technically this projection is accomplished by diagonalizing the unprojected Dirac 
Hamiltonian in a basis consisting of free Dirac spinors with momentum less than 
the cut-off A. That is, the Fock space for quarks on which the baryon is built, is 
truncated by the cut-off [T2[]. We will comment on the special role of the valence 
quark level later. The baryon number is defined as the asymmetry in the spectrum 
of h\. Apparently the just defined regularization yields unit baryon number for the 
functional (|i~2|). 

A self-consistent soliton solution would correspond to a local minimum of the 
energy functional (0). In contrast to the multi-quark systems discussed above we 
expect a strong coupling to isospin and hence the pions for a single baryon. Of 
course, this driving role of the pions is a well-known feature of the soliton picture 
for individual baryons where the pion fields assume a hedgehog form. We therefore 
choose the static ansdtze 

a{r) = a(r) and ir(r) = rix(r) (15) 

which introduces the two radial profile functions a(r) and 7r(r). These profile func- 
tions obey the self-consistency conditions 



N, 



air 
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An 



GjdQ {^/^va, - i £ tft/re \ (16) 



47T 



*(0 = -^Gldoj^z^r-f^-^^z^r-f^^ (17) 



with being the eigenfunctions of the projected Dirac Hamiltonian (|T3|). Suitably 
for the spherical configuration fll5|) the free basis to diagonalize h A is discretized by 
requiring appropriate boundary conditions for the quark wave-function at a finite 
but large radius, D. Of course, the existence of the continuum limit, D — ► oo is 
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crucial. States with momenta larger than A are omitted from that basis. In addition 
the basis states are labeled by the conserved grand spin quantum number with the 
valence level possessing grand spin zero. For more details on these technical aspects 



we refer to the literature [19|] 



The search for a self-consistent solution starts with choosing a trial configuration 
for the profile functions cr(r) and 7r(r). The resulting eigenstates of are subse- 
quently employed to update these profile functions according to the equations of 
motion (|l^,|l^). Iteration of this procedure eventually yields a stable configuration. 
Unfortunately we have been unsuccessful in finding a non-trivial configuration for 
the model parameters of table p] when taking the valence quark level to be the lowest 
eigenstate of the projected Hamiltonian. This includes the extreme case A = m, 
which supposedly yields strong binding (cf. table 1). This means that during the it- 
eration the meson profiles tend to the vacuum configuration a(r) = m and 7r(r) = 0. 
It should be noted that much care had to be taken in gaining this result as there 
appeared to be spurious finite size effects which faked a non-trivial solution. Al- 
though the typical extension of a soliton solution is expected to be of the order of 
lfm we were enforced to take D = 25fm or larger to avoid these finite size effects. 

On the other hand it should be noted that taking the valence level as resulting 
from diagonalizing the projected Hamiltonian h/y corresponds to also regularizing 
the valence level which is an un-common approach^ We have therefore also studied 
the case where the valence quark contribution to the equations of motion stems 
from diagonalizing the un-projected Hamiltonian. In that case something peculiar 
happens. Upon iteration the system tends to approach the un-stable (peaked) 
configuration observed in the proper-time scheme. This configuration in particular 
is characterized by a negative valence quark eigenenergy e va i. However, as e va i turns 
negative, it is considered part of the distorted vacuum and must be taken from the 
projected Hamiltonian, which, as discussed above, does not yield stable solutions. 
Hence in the iterative approach the configuration fluctuates about a point where e va i 
flips signF]. Apparently that is at best a saddle point solution of the equations of 
motion and cannot be considered a real solution. 



4 Renormalizable extension of the NJL model 

As already discussed in connection with the proper-time regularization scheme, the 
scalar field may develop an infinitely high and narrow peak at the origin, which 

4 In the proper-time scheme, for example, this would correspond to a non-integer baryon number 
and also destabilize the soliton. In the case of the three-momentum cut-off, however, demanding 
unit baryon number does not prohibit the regularization of the valence level. 

5 We have also considered the inclusion of a finite 'chemical potential' —m < fi < m against 
which we measure the single quark energies. In that case the configuration oscillates about the 
point e val = /i. 
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finally destroys the soliton. It is intuitively clear that any regularization scheme 
which does not suppress the high momentum components contained in that peak 
will encounter this problem. Here, we want to shed some light on how this instability 
occurs. The Wigner-Kirkwood expansion (section 2) proves suitable to give a simple 
explanation of the phenomenon. 

For definiteness we consider a non-trivial and renormalizable version of the NJL 
model, where the necessary counterterms, namely the mesonic kinetic energies and 
quartic mesonic interactions, are added to the Lagrangian ([!]). Here we give only a 
brief description of the renormalization procedure with particular emphasis on the 
soliton sector, for the details we refer to |16| |. For simplicity we consider the scalar 
field only, pseudo-scalar fields may be added straightforwardly 



£ = q{iYd. - cr)q + ^-d^Pa - 



1 2 Ao 
-a — 



a 



This amounts to a linear sigma model coupled to quarks which is known to be 



renormalizable ||17|| . According to the renormalization prescription JT6| the bare 
couplings indicated with subscripts zero are replaced by the renormalized (finite) 
parameters g a and The quadratic and logarithmic divergencies (denoted by 
I qua d and Ii g in [|T(J) stemming from the quark loop are isolated and expressed as 
momentum integrals, 



A = -iN c Ti log (i 7 M <9 M - a) + J d 4 xC M 



with 
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(2tt) 4 (p 2 - m 2 ) 3 
d A p (p 2 — 2m 2 ) 
(2tt) 4 {p 2 - m 2 ) 2 
d A p 1 



a 



(19) 



c 



(2vr) 4 (p 2 - m 



2\2 



a 



Ag 2 m? 

Note that the scalar field used here a = m + g a a' is connected to a' defined in 
and that the parameters are related /z 2 = 4m 2 (1 + Ncg 2 /Qii 2 ) such that the Nambu 
relation for the sigma mass m a = 2m holds. In principle the action (|T9|) is well 
suited to study the soliton in the renormalized extension of the NJL model. This is 
in particular the case because the counterterms, which are completely fixed in the 
meson sector, also render finite the fermion determinant in the soliton background. 
However, this is technically quite involved because identical regularization schemes 
have to be employed for the functional trace and the counterterms. In a numerical 
treatment this is technically quite complicatedF] and beyond the scope of the present 
paper. 



' Cf. ref EG] for a suitable path to approach this problem. 
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Instead we would like to gain some first insight by employing the much simpler 
semi-classical method described in section 2. It is noticed that the infinities of the 
extra terms contained in (^) just cancel those appearing in and (eq.§)) 
when A — ► oo (p*- -* and are finite in that limit). This is not obvious, apparently 
the WK expansion reproduces all singularities of the exact model correctly. Stated 
otherwise, the WK expansion allows us to carry out the renormalization program 
analytically and avoiding a numerical regularization of the fermion determinant. 

As the resulting action is a perfectly finite functional of the form (0) we can 
proceed analogously to subsection 3.1. We find that quark matter (symmetric phase 
with cr = 0) is in principle bound for g 2 > 24tt 2 /Nc- For finite systems the variation 
in the interior of the soliton leads to a non-linear differential equation similar to flTOl) . 
However, it turns out that in contrast to fllPP this differential equation possesses 
no stable solutions. Using a relaxation method for the solution of the non-linear 
differential equation with appropriate boundary conditions, the development of the 
infinitely high and narrow peak in the scalar profile at the origin is observed quite 
similar as reported in the proper-time regularization scheme |7], BJ. The kinetic part 
of the functional (@) inside the soliton (A = ep = Jp 2 F + a 2 ) as it follows from (|[), 



N, 



c 



1 

g 2 a 6tt 2 



- 3£n PF + eF ^ 



Pf 



m 



(W) 



(20) 



explains the cause of this result. Close to the center of the soliton, where pp increases 



towards ep 



m 



(er — > 0), the second term changes sign and eventually exceeds the 



constant and positive contribution from the first term. As a consequence the total 
kinetic energy becomes negative. Exactly this happened in all considered cases quite 
independently from the size of the soliton. Therefore, within the WK expansion, 
it is obviously the kinetic energy being no longer positive definite which causes the 
instability of the soliton. Of course, this is equivalent to the observation of an 
infinitely narrow and high peak in the scalar profile. 



5 Conclusions 

In this short note we have reported on the investigation of the soliton formation in 
the NJL model with local four quark interaction and various regularization proce- 
dures. Our interest was focused on solitons in the absence of additional constraints, 
which are commonly imposed in order to achieve stabilization. 

According to the investigations presented here the following picture emerges. 
Whether or not stable solitons may be found in the NJL model depends on the chosen 
regularization scheme. First of all there are the more sophisticated regularization 
schemes as e.g. Schwinger's proper-time regularization, which do not limit the high 
momentum components contained in the Dirac sea. Solitons in these schemes are 
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unstable because the scalar field develops an infinitely high and narrow peak at the 
origin. Using the WK expansion we were able to show that within a renormalizable 
extension of the NJL model the instability is related to the mesonic kinetic energy 
density not being positive definite. 

Apparently the three-momentum cut-off regularization scheme circumvents the 
above addressed problem by successfully suppressing the high momentum compo- 
nents contained in the sharp peak of the scalar field. This regularization scheme 
is therefore particularly appropriate to study the soliton formation in the uncon- 
strained NJL model. In detail our findings are: 

• Hadronic matter is bound for A/m < 1.80. Since finite systems add surface 
repulsion no solitons exist for larger ratios. 

• In order to describe finite size effects the lowest order gradient corrections 
to the Thomas-Fermi method in the WK expansion have been considered. 
We find that the binding energy per quark is continuously weakened with 
decreasing particle number until the soliton ceases to exist. This behavior, 
caused by the repulsive kinetic terms active in the surface of the soliton, reveals 
a peculiarity of the NJL model, which consequently prefers the formation of 
large clusters of hadronic matter. 

• The critical size, or equivalently the critical quark number of the soliton obvi- 
ously depends sensitively on the ratio A/m. We have computed the limiting 
function above which soliton formation takes place. Because of the semi- 
classical approximation employed, this result should not be trusted for small 
quark numbers. 

• An extensive search for a self-consistent configuration with unit baryon number 
was performed. No solitons were found, not even with an unphysically low 
cut-off A/m = 1 which should provide optimal attraction. In all cases the 
iterative procedure ran into the trivial configuration. Allowing the distinct 
valence quark not to undergo regularization, the sharp peak in the scalar field 
reappeared. 

Concludingly we have to state that although we do find solitons in the three- 
momentum regularized NJL model without further constraints, the results are un- 
satisfactory in various respects, in particular of course for the unit baryon number 
case. It seems that the NJL model together with an adopted regularization proce- 
dure alone cannot account for a reliable description of individual baryons as solitons. 
However, this may not necessarily be bad news for model builders. As discussed, 
the stabilizing extensions of the model can well be motivated by phenomenological 
considerations. 
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Finally, we would like to further add a few comments concerning the stability. 
Within the semi-classical approach the multi-quark solitons discussed in section 
3 are stable against variations of the scalar field. Independently we checked that 
the infinite system is at least locally stable also when spatially dependent pseudo- 
scalar fields are allowed. However, for finite systems, in particular of course for very 
small particle numbers, we may not exclude that in a self-consistent calculation 
pseudo-scalar fields leading to instabilities are excited. For the three-quark system 
we did not find a non-trivial solution. We may conjecture that the quartic meson 
self-interaction inherent in the NJL model is too weak in order to stabilize this 
soliton. 

On the other hand, with the chiral circle condition imposed, the proper-time 
regularized NJL soliton seems to be stable. Fluctuations in the lowest channels 



(grand spin) were searched and no instabilities detected [21]. It is natural to assume 
that this remains still valid when the chiral circle condition is softened and replaced 
by an additional quartic meson self- interaction. Thus, in both cases the hedgehog 
presumably represents the stable minimum configuration. 



For other versions of the model as e.g. its non-local extensions p2| |2^| the latter 



statement is not at all obvious. Although there the hedgehog solution is reported to 
be stable with respect to monopole deformations it is not excluded that this soliton 
leaks through other normal modes. 
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